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I. INTRODUCTION 

Single-inclusive hadron production is the most appro- 
priate and powerful benchmark for challenging altogether 
our understanding of the partonic structure of nucleons, 
the dynamics of hard Quantum Chromodynamics (QCD) 
interactions, the validity of QCD factorization, and the 
way in which quarks and gluons produce the detected 
final-state hadrons. The last few years have witnessed 
a remarkable improvement in both precision and vari- 
ety for this kind of measurements, which are expected 
to continue yielding crucial pieces of information in the 
future. 

In a recent article [l| , we have demonstrated the feasi- 
bility of performing a next-to-leading order (NLO) com- 
bined QCD analysis of single-inclusive pion and kaon 
production data, coming from electron-positron annihi- 
lation, semi-inclusive deep-inelastic lepton-nucleon scat- 
tering (SIDIS), and proton-proton collisions, collected 
over a wide kinematic range. The analysis not only al- 
lowed a very good global description of all these processes, 
but provided also accurate parametrizations for the 
parton-to-pion and parton-to-kaon fragmentation func- 
tions, which encode the details of the non-perturbative 
hadronization process relevant for the perturbative QCD 
(pQCD) framework based on the factorization theorem. 

In the following, we extend the analysis of Ref. ^|] first 
to the case of single- inclusive proton and anti-proton pro- 
duction, again including not only the electron-positron 
annihilation data [2, iJ, 3, [E, @ used in all previous fits 
0)11, Hi J but, for the first time, also the very precise data 
recently obtained in proton-proton collisions at RHIC 
p^ . where final state protons and anti-protons are dis- 
criminated. 

Next, after having obtained reliable fragmentation 
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functions for the three lightest and most copiously pro- 
duced charged hadron species, pions and kaons from 
Ref. [l| and (anti-)protons, we analyze inclusive, i.e., 
unidentified, charged hadron yields obtained in electron- 
positron annihilation [2 J3, H, [ sl [III [H, [H, [l3j_proton- 
(anti-)proton collisions[lE [ll,ll37and SIDIsfls'l. Here 
we are aiming at an extraction of the contribution from 
the "residual" charged hadrons other than pions, kaons, 
and (anti-)protons, to the inclusive charged hadron frag- 
mentation functions. 

Only from a global QCD analysis we can obtain in- 
dividual quark and anti-quark fragmentation functions 
for unidentified positively or negatively charged hadrons 
without having to make assumptions on the relation be- 
tween favored (valence-type) and unfavored (sea-type) 
contributions. Such assumptions have been shown to 
be often not adequate for reproducing presently avail- 
able single-inclusive hadron production yields beyond 
those obtained in electron-positron annihilation In 
addition, data from proton-(anti-)proton collisions are 
crucial for reducing the uncertainties in the gluon-to- 
hadron fragmentation function 1] since scaling violations 
of e+e~ data are not much of a constraint due to the lack 
of precision data at more than one value of the center-of- 
mass system (c.m.s.) energy ^/s. Both resulting new sets 
of NLO fragmentation functions reproduce with remark- 
able accuracy the data available and hence complement 
consistently our previous studies for pions and kaons. 

Previously available sets of inclusive charged hadron 
fragmentation functions are all based on fits to e+e~ data 
only and either sum up the results obtained for pions, 
kaons, and (anti-)protons, ignoring possible contributions 
from heavier charged hadrons ?,'8] or the residual hadron 
contribution already includes protons and anti-protons 
(Tot . The analysis of Ref. provides only fragmen- 
tation functions for unidentified charged hadrons, again 
based on e"'"e~ data only, however, an attempt was made 
to further constrain the gluon fragmentation function by 
comparing to one of the many data sets from proton- 
(anti-)proton collisions. 
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In the next section, we briefly highlight the main fea- 
tures of the present analysis, specifically discussing our 
choice of parametrizations for the proton and the residual 
charged hadron fragmentation functions and emphasiz- 
ing the charge conjugation and flavor symmetry assump- 
tions we still have to make. In Section III we present our 
results for both NLO fits, and assess the uncertainties 
involved with the help of the Lagrange multiplier tech- 
nique. In Sec. IV we summarize our main results. For 
completeness, the Appendix contains the results of the 
global analyses performed at leading order (LO) accu- 
racy, which compares significantly less favorable to data 
than our NLO fits. 



II. 



OUTLINE 



In the present analysis we work in the well established 
framework of pQCD for single-inclusive hadron produc- 
tion processes at NLO accuracy, thoroughly discussed 
and implemented in Ref. We make an extensive use 
of the Mellin transform technique 2l| developed for a 
fast computation of the exact NLO cross sections in each 
step of a global minimization procedure. We refer the 
reader to |l| and references therein for further technical 
details on the general QCD framework and the assess- 
ment of uncertainties in global fits with the help of La- 
grange multipliers. The analysis proceeds in two stages, 
the first dedicated to the extraction of the parton-to- 
(anti-)proton fragmentation functions, followed by the 
determination of inclusive charged hadron fragmentation 
functions in the second. 



A. Proton fragmentation functions 

In the first step of our global analysis we aim to deter- 
mine individual fragmentation functions Df^ for quarks 
and anti-quarks of all flavors i as well as for gluons into 
either protons {H — p) or anti-protons {H = p). At vari- 
ance with previous analyses the use of data that 
discriminate between p and p in the final-state, allows 
us to extract fragmentation functions for either particle. 
As explained in Ref. in order to have the flexibility 
required by charge separated distributions and to accom- 
modate the additional data, we adopt a somewhat more 
versatile functional form for the input distributions at 
scale ^0 than in 0, H, Q 



iV,z"'(l-^)ft[l + 7,(l-z)^ 
B[2 + ai, I3i + l]+ 7,B[2 + a„ 13, 4 



r,(l) 



where B[a^h] represents the Euler beta function and Ni 
is normalized such to represent the contribution of to 
the second moment dzzDf^ {z, /zq) entering the energy- 
momentum sum rule, z denotes the fraction of the frag- 
menting parton's energy taken by the produced hadron 



H. As in [H, the initial scale /io in Eq. ^ for the Q^- 
evolution is taken to be /io = 1 GcV for the light u, d, s 
partons and the quark masses for the heavier ones. 

The parameters describing the fragmentation functions 
D^{z,fio) in Eq. ([1]) are determined by a standard 
minimization for K data points, where 



K 



(2) 



Ej is the measured value of a given observable, 5Ej the 
error associated with this measurement, and Tj is the 
corresponding theoretical estimate for a given set of pa- 
rameters in ([1]). Since the full error correlation matrices 
are not available for most of the data entering the global 
analysis, we take the statistical and systematical errors 
in quadrature in 5Ej. 

In order to reduce the number of parameters in ([1]) to 
those that can be effectively constrained by the data, we 
impose, as usual, certain relations among the individual 
fragmentation functions. We have checked in each case 
that relaxing these assumptions indeed does not signifi- 
cantly improve the total of the fit ^ to warrant any 
additional parameters. In detail, for the fragmentation 
of light quark fiavors to a proton we assume the same 
shape for up and down quarks and for up and down anti- 
quarks with the same ^-independent normalization ratios 
N, i.e., 



and D 



p _ 



(3) 



The relation between quark and anti-quark fragmenta- 
tion functions for u and d flavors is determined by the 
global fit through 



(4) 



with f3 constrained to be positive. For strange quarks it 
is assumed that 



(5) 



with the SU(3)-breaking parameter N' independent of z. 
For charm, bottom, and gluon-to-proton fragmentation 
we find no improvement in the total for 7i 7^ in ([1]) , 
hence we set 7^ = 0. 

To obtain the corresponding fragmentation functions 
for anti-protons we assume charge conjugation sym- 
metry, i.e.. 



(6) 



leaving in total 17 free parameters to be determined by 
the global fit. 



B. Charged hadron fragmentation functions 

The fragmentation functions for unidentified pos- 
itively charged hadrons {H = h^) are defined by 



D 



DP + 



(7) 
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where and Df^ were already determined in Rcf. 1] . 
jjres+ denotes the residual contribution of positively 
charged hadrons other than pions, kaons, and protons to 
the inclusive sum . A definition analogous to Eq. ([7]) 
is used for the fragmentation into negatively charged 
hadrons. 

Since pions are much more copiously produced than 
heavier kaons and protons are even less abundant, it is 
natural to expect that Eq. ([7]) is strongly dominated by 
of +D^^ +D^^ leaving Df^*"^ to be small. Nevertheless 
it is an important consistency check to extract a small but 
non-vanishing Dl^^ (and Dl"^^ ) from data to actually 
confirm this hierarchy. If the global fit would require, for 
instance, a large or even a negative residual contribution, 
the usefulness of the previously extracted fragmentation 
functions for light hadrons would be in jeopardy. 

For Dl'^'^^ it turns out that the data are most econom- 
ically described by assuming full SU(3) fiavor symmetry 
for both quarks and anti-quarks 

and 

Dr^ = Dr^ = Dr\ (9) 

respectively, in the ansatz ([1]). Again, 

2D-^+ (10) 

and for D^''^^ we set jc,b = in ([1]), however, has 

a preference for ^ 0. For we assume charge 

conjugation symmetry, i.e., 

leaving in total 18 free parameters to be fitted here. 

We note that the Q^-dependence for all fragmentation 
functions is computed with the appropriate NLO evolu- 
tion equations as explained in detail in Ref. [l|. Uncer- 
tainties in the extraction of the parameters in Eq. ([T]) 
will be assessed with the help of the Lagrange multiplier 
technique as described again in |l| and references therein. 

III. RESULTS OF THE GLOBAL ANALYSES 

In this section we discuss in detail the results of our 
global analyses of fragmentation functions [3l| for pro- 
tons and residual charged hadrons as outlined above. We 
present the parameters of the optimum fits describing 

the NLO Df'^^" at the input scale no and compare to 
all data sets used in the analyses, including individual 
values. Detailed comparisons are made with previous fits 
based exclusively on the available electron-positron an- 
nihilation data in Refs. 3, @, and [l^, in the following 
labeled as AKK, HKNS, and KRE, respectively. 



A. NLO analysis of proton fragmentation functions 

Three different processes allow the extraction of the 
(anti-)proton fragmentation functions in our global anal- 
ysis. First of all, we have the "standard" electron- 
positron annihilation data @j B S Q customarily in- 
cluded also in most of the previous analyses 0, H, • A 
characteristic feature of e+e" annihilation data in gen- 
eral is that they only provide information on a certain 
hadron species summed over the charges, i.e., in this 
case on the sum of protons and anti-protons. Besides 
the fully inclusive measurements, SLD Q and DELPHI 
[B| also give results for "fiavor enriched" cross sections, 
distinguishing between the sum of light u, d, s quarks, 
charm, and bottom events. The quark flavor is deter- 
mined from Monte-Carlo simulations and refers to the 
primary qq pair created by the intermediate Z-boson or 
photon. Such results have to be taken with a grain of salt 
as they cannot be unambiguously interpreted and calcu- 
lated in pQCD. Nevertheless, in our analysis of pion and 
kaon data Q we always found good agreement between 
data and theory. 

In Fig. [1] we provide a detailed comparison in terms of 
"(data-theory) /theory" between all the e+e" data sets 
included in the analysis and our NLO results. The best 
agreement is found for ALEPH and DELPHI inclusive 
measurements, although SLD and TPC, as well as fiavor 
tagged data sets, also agree well with the fit within the 
fairly large experimental uncertainties. The fiuctuations 
in the data are significantly larger than in the case of 
pions, but roughly comparable to those found for kaons, 
see Ref. ^1]. Other recent fits [1, all reproduce the 
data shown in Fig.[l]equally well within the experimental 
uncertainties. Since the current e'^e" data basically only 
constrain the total singlet fragmentation function 
at scale Mz, all sets differ considerably in their individual 
flavor content as we shall see below. 

We wish to recall that the range of applicability for 
fragmentation functions is severely limited to medium- 
to-large values of z, see, e.g., Ref. [l[. In order to avoid 
the potentially problematic low-^ region, we exclude from 
the fit all data with energy fractions lower than Zmin — 
0.1. The extrapolation of our fit, however, reproduces 
the trend of the data reasonably well also below Zmin as 
indicated in Fig. [T] 

The second key ingredient in our global analysis are 
the single-inclusive hadron production data from proton- 
proton colhsions at BNL-RHIC taken by STAR ^ at 
mid-rapidity \y\ < 0.5 and shown in Fig. [21 These data 
discriminate final-state protons from anti-protons and 
hence, in principle, allow to separate quark-to-proton 
and anti-quark-to-proton fragmentation functions in the 
fit. However, at the presently accessible range of trans- 
verse momenta pt and at mid-rapidities the production 
of single-inclusive hadrons is mainly driven by gluon- 
induced processes and fragmentation, turning these data 
into the best constraint on the gluon fragmentation func- 
tion DP at large values of z currently available. Figure 
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FIG. 1: "(data-theory) /theory" comparison of our NLO results for the electron-positron annihilation cross section into protons 
and anti-protons with the data sets used in the fit, see also Tab. ITT] 
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FIG. 2: Upper panels: comparison of our NLO results (solid lines) for single-inclusive proton (anti-proton) production pp — » 
p{p)X at y^s = 200 GeV with STAR data [l^I at mid-rapidity \y\ < 0.5 using fj,f = fj,r = pr- Also shown are the results 
obtained with the HKNS 9^ parametrization. The shaded bands indicate theoretical uncertainties when all scales are varied in 
the range pt/2 < fJ^f ~ fJr < 2pT. Lower panels: "(data-theory) /theory" for our NLO results. 
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FIG. 3: Comparison of the OPAL "tagging probabilities" 0] 
for protons plus anti-protons, as a function of the minimum 
Xp, with our NLO results (solid lines). Also shown are the 
results obtained with the AKK parameterization (dashed 
lines) . 

[5] shows the results of the fit compared to the data. The 
theoretical uncertainties related to a particular choice of 
the factorization and renormalization scales, ^/ and fir, 
respectively, in the computation of the NLO cross sec- 
tion, indicated by the shaded bands, are non- negligible. 
This has been taken into account in the minimiza- 
tion procedure as a conservative additional 5% relative 
uncertainty. 

We have checked that by excluding data at low pT 
values, Pt ^ 3 GeV, the results of the global fit remain 
essentially unchanged. The fact that the pt spectrum is 
well reproduced down to 1 GeV, the mass of the proton, 
may turn out to be accidental, but this can only be veri- 
fied once data at higher px become available to map out 
the Pt slope more precisely. We note that the BRAHMS 
data [22] at very forward rapidities y ~ 3 and pt ^ 4 GeV 
have not been included in our global analysis. These data 
show a very pronounced difference between the p and p 
yields which cannot be understood. The origin of the 
excess of protons over anti-protons by a factor of about 
10, not observed by STAR at less extreme kinematics, 
remains an open question [221] . 

Finally, a third source of information on the parton-to- 
proton fragmentation is provided by the OPAL "tagging 
probabilities" rjf [6j, also included in the AKK fit [8] 
but not in 9] , which are sensitive to the fiavor of q + q 
fragmentation functions. The rif{xp) represent "proba- 
bilities" for a quark flavor i to produce a "jet" contain- 
ing the (anti-)proton with a momentum fraction z larger 
than Xp. Figure [3] compares the results from our and 
the AKK fit with the data. As discussed at length in 
Ref. [ij and also mentioned in the discussion of the fiavor 
tagged SLD and DELPHI data above, any flavor tagged 
information is highly model dependent and difficult to 
interpret within pQCD beyond the LO approximation 
involve non trivial theoretical uncertainties. To take this 
into account, we assign an up to 10% extra normalization 



uncertainty to the rjf data in the minimization. The 
agreement between OPAL data and theory is reasonably 
good. 

Tables U and |TT1 show the values obtained for the pa- 
rameters in Eq. ([1]) specifying the optimum flt of proton 
fragmentation functions D^(z,/io) at NLO accuracy and 
summarize the values for each individual set of data 
included in the global analysis, respectively. In cases 
where the normalization uncertainty of the experiment 
is not included in the error bars of the data, we apply 
a free normalization factor constrained to vary within 
range quoted by the experiment. The values of the nor- 
malization factors resulting from the flt are also included 
in Tab. HH 
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FIG. 4: Upper panels: individual fragmentation functions for 
protons zDf{z, Q^) at = M%. Middle panels: ratios of our 
fragmentation functions to the ones of HKNS [9| . The dotted 
line indicates the ratio for singlet combination of fragmenta- 
tion functions 2; _D|.. Lower panels: ratios of our fragmentation 
functions to the ones of AKK [|. 

Figure [4] shows the resulting new set of fragmentation 
functions for protons for different flavors, evolved to the 
scale of the mass of the Z-boson, and in the middle and 
lower panels, comparisons to the fits of HKNS [oj and 
AKK Q , respectively. As expected, the present fit agrees 
well in the singlet combination with previous extrac- 
tions based only on electron-positron annihilation data, 
but we find significant differences in the gluon fragmen- 
tation function, constrained in our fit not only by the 
scale dependence of the data but also by the STAR data 
shown in Fig. [21 In addition, there is a sizable difference 
between our and that of HKNS which is mainly a con- 
sequence of the SU(3) flavor symmetry imposed for the 
unfavored fragmentation functions in the latter analysis 

a- _ 

As in Ref. [l| , we make use of the Lagrange multiplier 
technique [2^ in order to give a representative picture of 
the typical uncertainties characteristic of the fragmenta- 
tion functions obtained from the global flt. In Fig. [H] we 
show the x^-profiles as a function of the range of varia- 
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TABLE I: Parameters describing the NLO fragmentation 
functions for protons, Df{z,^o), in Eq. ^ at tfie input scale 
/io = 1 GeV. Inputs for the charm and bottom fragmenta- 
tion functions refer to = mc — 1.43 GeV and /io — mt — 
4.3 GeV, respectively. 



flavor i 






ft 


li 


5r 


u + u 


0.091 


-0.222 


1.414 


15.0 


3.29 


d + d 


0.058 


-0.222 


1.414 


15.0 


3.29 


u 


0.034 


-0.222 


2.024 


15.0 


3.29 


d 


0.022 


-0.222 


2.024 


15.0 


3.29 


s + s 


0.043 


-0.222 


2.024 


15.0 


3.29 


c + c 


0.076 


-0.899 


5.920 


0.0 


0.00 


b + b 


0.044 


-0.034 


10.000 


0.0 


0.00 


9 


0.014 


6.000 


1.200 


0.0 


0.00 



TABLE II: Data used in the NLO global analysis of proton 
fragmentation functions, the individual values for each set, 
the fitted normalizations, and the total of the fit. 



experiment 


data 


rel. norm. 


data points 






type 


in fit 


fitted 




TPC [2] 


inch 


1.06 


12 


7.7 


SLD [3] 


inch 


0.983 


18 


12.0 




"lids tag" 


0.983 


9 


10.5 




"c tag" 


0.983 


9 


9.8 




"& tag" 


0.983 


9 


8.9 


ALEPH [4| 


inch 


0.97 


13 


11.5 


DELPHI [5] 


inch 


1.0 


12 


3.9 




"wds tag" 


1.0 


12 


0.6 




"6 tag" 


1.0 


12 


9.1 


OPAL [6] 


"m tag" 


1.10 


5 


7.6 




"d tag" 


1.10 


5 


13.5 




"s tag" 


1.10 


5 


5.0 




"c tag" 


1.10 


5 


4.9 




"fe tag" 


1.10 


5 


5.5 


STAR [10] 


P 


0.95 


14 


35.4 




P 


0.95 


14 


26.0 


TOTAL: 






159 


171.9 



tion of the truncated second moments of the individual 
fragmentation functions of flavor i, 

i^P{xp,Q')= [' zDP{z,Q')dz, (12) 

J Xp 

for Xp — 0.2 and Q — b GeV, around the values obtained 
for them in the best fit to data, Tyfp. Roughly speak- 
ing, the typical uncertainties in the second moments of 
the proton fragmentation functions range between 20% 
and 25%, allowing for a conservative increase Ax^ of 2% 
in the total of the fit, except for the gluon and bot- 
tom fragmentation functions, where the uncertainties are 
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FIG. 5: Profiles of x^ for the NLO proton fragmentation 
fit as a function of the truncated second moments r]f{xp = 
0.2, Q'^ — 25GeV'^) for different flavors. The moments are 
normalized to the value rjfg they take in the best fit to data. 

closer to 10%. The rather stringent constraint for bot- 
tom comes from the availability of DELPHI and SLD 
flavor tagged data, while for charm only SLD data with 
larger errors are available. Note that neither charm nor 
bottom fragmentation play any role in the description of 
the STAR data due to the relatively low scales of 0{pt) 
involved. For the gluon fragmentation function D^, the 
STAR data make the difference. We wish to stress that 
compared to pions and kaons the parton-to-proton 
fragmentation functions are much less constrained at 
the moment. This is also reflected in the stronger as- 
sumptions, Eqs. ©-(in]), which have to be imposed on 
the fit in order to be able to determine all parameters. In 
particular the lack of (anti-)proton production data from 
SIDIS prevent a more reliable separation of favored and 
unfavored fragmentation functions. Also the gluon frag- 
mentation is currently mainly determined from a single 
set of data (STAR). Here, possible future high precision 
data from B-factories would open up the possibility for 
studies of scaling violations in e+e^ annihilation, which 
should help to further constrain D^. 



B. NLO analysis of charged hadrons fragmentation 
functions 

As outlined in Sec. IIIBl rather than extracting frag- 
mentation functions for inclusive charged hadrons from 
scratch, we take advantage of the knowledge already ac- 
quired on pion, kaon, and proton fragmentation func- 
tions. In the following, we isolate and determine the 
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FIG. 6: Left-hand side: comparison of our NLO results for the electron-positron annihilation cross section into inclusive charged 
hadrons with the data sets used in the fit, see Tab. II VI Right-hand side: (data-theory) /theory for our NLO results for each of 
the data sets. 




FIG. 7: Same as in Fig. [S] but now for light ("uds") and heavy quark ("c,b") tagged cross sections. 
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contributions coming just from residual charged hadrons 
Dl'"' as defined in Eq. In this way, there is no 

need to make assumptions on the relations between the 
inclusive charged hadron fragmentation functions - nec- 
essary for fitting purposes - which may hold for one of 
the most abundant hadron species like pions, but not nec- 
essarily for the others. Factoring out the contributions 
from charged pion, kaon, and (anti-)proton fragmenta- 
tion functions, we are left with a comparatively small 
residual contribution Dp** , whose flavor symmetry as- 
sumptions, see Sec. Ill Bl have very little impact on the 
inclusive charged hadron fragmentation functions and the 
quality and reliability of the fit. 

We begin by discussing the remarkable measurements 
of single-inclusive charged hadron production in e^e~ an- 
nihilation, collected by a variety of experiments 0, d, 0, 
SUiljEI- Figures [5] and [7| show detailed comparisons be- 
tween the available data and the results of adding the al- 
ready known pion, kaon, and proton fragmentation func- 
tions to the outcome of the fit for the residual hadron 
contributions Dl"^^ . 

Tables Hill and HVl show the values obtained for the pa- 
rameters in Eq. ([1]) specifying the optimum fit of residual 
charged hadron fragmentation functions Dl"^"^ {z , ^q) at 
NLO accuracy and summarize the values for each in- 
dividual set of data included in the fit, respectively. As 
in Tab. |TT1 for each set of data, values of the normal- 
ization factors resulting from the fit and constrained to 
vary within range quoted by experiment, are also given 
in Tab. |TVl 

All data sets shown in Figs. [5] and [7] are very well re- 
produced by our new fit, again closely following the trend 
of the data again also below Zmin = 0.1. This is most 
apparent in the (data-theory) /theory comparisons. The 
impressive quality of the fit serves as an important cross- 
check of the consistency of entire fitting procedure for all 
hadron species since the main contributions to the inclu- 
sive charged hadron cross sections in Figs. [H] and [7] are 
already fixed by the pion, kaon, and proton fits. For in- 
stance, overshooting the data by just summing up pion, 
kaons, and protons would have pointed to some serious 
inconsistency in the global analysis. As expected, we also 
find that the residual contribution from D^*^* is indeed 
small and - another non-trivial check - positive, see be- 
low. Note that we estimate an average uncertainty of 
5% in all theoretical calculations of the inclusive charged 
hadron cross sections stemming from the propagation of 
uncertainties of pion, kaon, and proton fragmentation 
functions. This additional uncertainty is included in the 
minimization procedure for determining D^'^'^ . 

As can be inferred from Tab. IIVI our global analysis 
contains also several other data sets which help to further 
constrain certain aspects of the residual charged hadron 
fragmentation functions and, more importantly, provide 
additional consistency checks. Longitudinal cross section 
measurements in e^e~ annihilation are expected to help 
to constrain giuon fragmentation since, contrary to in- 



TABLE III: Parameters describing the NLO fragmen- 
tation functions for positively charged residual hadrons, 
Dl'^"'^ (z, fio), at the input scale fj,o = 1 GeV. Inputs for 
the charm and bottom fragmentation functions refer to fio = 
rric — 1.43 GeV and ~ mb ~ 4.3 GeV, respectively. 



flavor i 






ft 


7» 




u + u 


0.0038 


10.000 


1.20 


0.0003 


18.51 


u 


0.0001 


10.000 


21.20 


0.0003 


18.51 


c+c 


0.0752 


0.406 


3.91 


0.0000 


0.00 


b + b 


0.0936 


-0.150 


3.61 


0.0000 


0.00 


g 


0.0001 


-0.497 


9.99 


20.000 


14.75 



TABLE IV: Data used in the NLO global analysis of residual 
charged hadron fragmentation functions, the individual 
values for each set, the fitted normalizations, and the total x^ 
of the fit. 



experiment 


data 


rel. norm. 


data points 






type 


in fit 


fitted 




TPC [2] 


inch 


1.027 


17 


1.7 


SLD [3] 


inch 


1.006 


21 


13.2 


ALEPH [4] 


inch 


1.027 


27 


27.0 


DELPHI [5] 


inch 


1.0 


12 


6.8 




"uds tag" 


1.0 


12 


7.6 




"& tag" 


1.0 


12 


4.9 


TASSO [11] 


inch (44 GeV) 


1.0 


14 


11.5 




inch (35 GeV) 


1.0 


14 


19.9 


OPAL [12] 


inch 


1.0 


12 


5.7 




"uds tag" 


1.0 


12 


11.9 




"c tag" 


1.0 


12 


7.4 




"b tag" 


1.0 


12 


3.9 


ALEPH 14] 


long, inch 


1.0 


11 


1.8 


OPAL 


long, inch 


1.0 


12 


3.3 


DELPHI [14] 


long, inch 


1.0 


12 


11.6 




long, "uds tag" 


1.0 


12 


35.1 




long. "6 tag" 


1.0 


12 


4.7 


EMC [18] 


h+ 


0.987 


98 


99.1 




h- 


0.987 


99 


156.8 


CDF 115] 


630 GeV 


1.1 


16 


103.4 




1.8 TeV 


1.1 


37 


112.7 


UAl \16\ 


200 GeV 


1.1 


31 


111.5 




500 GeV 


1.1 


32 


44.5 




630 GeV 


1.1 


41 


214.2 




900 GeV 


1.1 


44 


118.1 


UA2 117] 


540 GeV 


1.1 


27 


89.3 


TOTAL: 






661 


1227.6 



elusive total cross sections, a non- vanishing 0{as) gluon 
coefficient function is already present at LO. For our NLO 
fits of the longitudinal cross section, we include the 0{al) 
coefficient functions ^]. In Fig. [8]we compare the result 
from the NLO global analysis for the longitudinal inclu- 
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FIG. 8; Same as in Figs. |6] and [7] but now for the longitudinal cross sections. 
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FIG. 9: Comparison of our NLO results (solid lines) for the positively and negatively charged hadron multiplicities in SIDIS, 
(l/A'^Dis)dA'''' /dzdQ^, with data from EMC [3 for different bins in z. Upper, middle and lower panels correspond to beam 
energies of 120, 200, and 280 GeV, respectively. Also shown are the results obtained with the KRE [l^ parameterization 
(dashed lines). 
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FIG. 10: Upper panels: comparison of our NLO results for single-inclusive charged hadron production pp — > hX, where 
h — (ft^ + h~)/2, with data from UAl [l^ and UA2 [l7|] for various c.m.s. energies y^s using fj,f = fir ~ Pt- The UAl data 
cover jyj < 2.5, except for 630 GeV, where \y\ < 3.0 and also h — (/i^ + h~). For UA2 the rapidity range 1.0 < |y| < 1.8 is 
covered. Also shown are the results obtained with the AKK [3| parametrization including only pious, kaons, and (anti-)protons. 
The shaded bands indicate theoretical uncertainties when all scales are varied in the range pt/2 < Hf = fir < 2pT. Lower 
panels: "(data-theory) /theory" for our NLO results. 



sive and flavor tagged cross section with the available 
LEP data 0, [H, [ij . With the exception of the ALEPH 
data, the statistical precision is not really impressive, in 
particular, in the fitted region z > Zmin, but all data are 
reasonably well described by the fit. 

Data on charged pion and kaon multiplicities in SIDIS 
have provided valuable information on scale dependence 
and the charge and flavor separation of the corresponding 
fragmentation functions in our previous analysis [l| . For 
inclusive charged hadrons, the data collected by EMC 
[l8| cover a rather large kinematic region with very nar- 
row (5^-bins, what makes it particularly suitable for a 
QCD fit. As expected, these data prove to be very valu- 
able in our global analysis here. Figure [S] shows the out- 
come of the fits compared to the EMC multiplicities for 
positively and negatively charged hadrons covering all 
beam energies of 120, 200, and 280 GeV. The good pre- 
cision of the data over the entire z range helps to further 
constrain the high z behavior of the quark and anti-quark 
fragmentation functions. An interesting thing to notice 
is that even though the KRE sets for pions and kaons 
overestimate the corresponding SIDIS multiplicities as 
demonstrated in [l[ , this effect is somehow compensated 
for in the charged hadron multiplicities shown in Fig. [5] 
at the expense of comparatively small residual fragmen- 
tation functions. In order to account for the error intro- 
duced by the finite size of the bins in z and Q^, as well 



as the uncertainties introduced by the use of the pion, 
kaon, and proton fragmentation functions, we include an 
additional 10% theoretical uncertainty in the mini- 
mization. There is also an additional 11% experimental 
normalization uncertainty not included in the errors 
bars shown in Fig. [Sj 

The final ingredient to our global analysis of fragmen- 
tation functions for unidentified charged hadrons are 
single-inclusive hadroproduction data from pp collisions 
measured by the UAl [l^ and UA2 [l3| collaborations at 
CERN, and by the CDF collaboration [H] at Fermilab's 
TeVatron. We have not included proton-proton collision 
data from fixed target experiments, since the validity of 
fixed order calculations in pQCD at lower energies is seri- 
ously in doubt [25]. The data span a range of c.m.s. ener- 
gies ^/i from 200 GeV to 1.8 TeV but do not discriminate 
different hadron charges. Due to the dominance of gluon- 
induced processes at the available small-to-medium val- 
ues of the hadron'spT, they mainly probe the gluon frag- 
mentation function . To give an example, in the case 
of UAl at = 630 GeV, where the available data cover 
the largest range of transverse momenta pt from 1 GeV 
to 22 GeV, in which the gluonic contribution decreases 
from around 90% in the lowest pr bin to a sizable 40% 
at the highest transverse momentum. Figures [10] and 1111 
show the comparisons between the results of our NLO fits 
to the single-inclusive cross sections and the experimen- 
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FIG. 11: Similar to Fig. [TO] but now comparing our NLO 
results to the CDF data at two different c.m.s. energies ^/s. 
The data are for/i = {h+ + h-)/2and cover \y\ < 1 in rapidity. 



tal data from CERN and Fermilab, respectively. Within 
the fairly large theoretical scale ambiguities, indicated 
by the shaded bands, which refer to varying the factor- 
ization and renormalization scales simultaneously in the 
range Pt/2 < /i/ = /ir ^ 2pT, the overall agreement be- 
tween theory and data is reassuringly good. The obtained 
gluon fragmentation functions into charged pions, kaons, 
and (anti-)protons, which dominate the single-inclusive 
cross sections shown in Figs. [TU] and [TTJ are consistent 
with data, and the residual contribution Dg'^'^ is very 
small. Again, we include a 5% theoretical error due to 
the choice of the scale in the minimization and another 
5% associated with the propagation of the uncertainties 
of pion, kaon, and proton fragmentation functions. 

The resulting fragmentation functions for positively 
charged hadrons , evolved to the scale Mz, are shown 
in Fig. [T^ and compared to the distributions of KRE [l^ 

and AKK "8]. Recall that in the AKK analysis £>f^ 
is approximated by the sum of pion, kaon, and proton 
fragmentation functions. Compared to sets of KRE and 
AKK, we find again good agreement in the total singlet 
contribution D'^ , but some noticeable differences in the 
charge and flavor separation, for instance, in D^^ , and in 
the gluon fragmentation D^^ . In general, the differences 
become larger towards z — > 1, which has been already ob- 
served for the pion, kaon [l[, and proton fragmentation 
ftmctions (Fig. |4|) . 

One striking feature of the distributions shown in 
Fig. [T^] is the only slightly broken flavor democracy for 
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FIG. 12; Upper panels; individual fragmentation functions 
for positively charged hadrons zD^ {z,Q^) at — M| for 
i — u + u, 2u, s + s, g, c + c, and 6-1-6. Middle panels: ra- 
tios of our fragmentation functions to the ones of KRE fl9| . 
The dotted line indicates the ratio for singlet combination of 
fragmentation functions zD^ . Lower panels: ratios of our 
fragmentation functions to the ones of AKK [^; note that 
is not available in the AKK analysis. 
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FIG. 13; Partial contributions of pions, kaons, and protons 
to the charged hadron fragmentation functions as a function 
oiz at = M|. 



q + q fragmentation functions into charged hadrons. To 
understand this better, the partial contributions of pions, 
kaons, and protons to the charged hadron fragmentation 
functions are shown in Fig. 1131 As can be noticed, the u 
quark and gluon fragmentation functions are completely 
dominated by the pion contribution, except for very large 
momentum fractions z, where proton-proton scattering 
data require some increase in the kaon and proton part. 
The strange quark fragmentation function for charged 
hadrons in the valence region is, as expected, dominated 
by kaons, while both kaons and pions contribute with a 
similar amount to the charm (and bottom) distributions. 
The residual contribution becomes sizable only for the 
heavy quark fragmentation functions. 

Even though the relative uncertainties of the residual 
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0.99 1 1.01 0.99 1 1.01 




0.98 1 1.02 0.98 1 1.02 




FIG. 14: Profiles of for the NLO residual charged hadron 
fragmentation fit as a function of the truncated second mo- 
ments T]T'^ [x^ = 0.2, = 50GeV^) for different flavors. 
The moments are normalized to the value r^io*^ they take 
in the best fit to data. The upper horizontal scales show 
variations relative to the change in the moments of the total 
charged hadron fragmentation functions . 

charged hadron fragmentation functions D\'^^^ are rather 
large, exceeding even 200% within a conservative 2% in- 
crease Ax^ in the total of the fit, the effect on i'f ^ is 
fairly small and amounts to less than about 2%. This is 
readily explained by the relative smallness of the residual 
charged hadron contribution to the sum in Eq. ([7]). The 
best constrained residual densities are those for charm 
and bottom, which are also the most significant ones in 
absolute size, although almost innocuous for SIDIS and 
pp data. In Fig. [T3]we show the x^-profiles of the global 
fit as a function of the range of variation of the trun- 
cated second moments of the residual charged hadron 
fragmentation functions Z3['^*^, see Eq. (fT^ . In the up- 
per horizontal scales we also show the impact of those 
variations relative to the change in the moments of the 
total charged hadron fragmentation functions . 



IV. CONCLUSIONS 

By extending our previous global analyses for pions 
and kaons to the case of (anti-)protons and unidenti- 
fied charged hadrons, we have completed a comprehen- 
sive study of single-inclusive hadron production within 
pQCD at NLO accuracy. Specifically, we have demon- 
strated that the pQCD framework based on factorized 
cross sections can consistently account for a large variety 



TABLE V: As in Tab. |T] but now describing the LO fragmen- 
tation functions into protons, D^{z,fj,o), at the input scale. 



flavor i 


Ni 






Si 


u + u 


0.094 


0.041 1.485 


15.0 


3.44 


d + d 


0.059 


0.041 1.485 


15.0 


3.44 


u 


0.036 


0.041 1.998 


15.0 


3.44 


d 


0.022 


0.041 1.998 


15.0 


3.44 


s + s 


0.045 


0.041 1.998 


15.0 


3.44 


c -I- c 


0.079 


-0.887 5.436 


0.0 


0.00 


h + b 


0.047 


0.103 10.00 


0.0 


0.00 


9 


0.029 


6.000 1.200 


0.0 


0.000 


TABLE VI: 


As in Tab.l 


TI]but now for the LO fragmentation 


functions into positively charged residual hadrons D^"^ 


+ 


flavor i 




Qi ft 


7» 




u + u 


0.0002 


9.986 1.543 


15.0 


19.90 


u 


0.0001 


9.986 21.543 


15.0 


19.90 


c+c 


0.147 


-0.051 2.792 


0.0 


0.00 


h + b 


0.113 


-0.574 2.949 


0.0 


0.00 


9 


0.0001 


-0.499 10.000 


20.0 


19.64 


TABLE VII: Same as 


in Tab, im but now 


at LO accuracy. 


experiment 


data 


rel. norm. data points 


2 

X 




type 


in fit 


fitted 




TPC [2] 


inch 


1.043 


12 


7.5 


SLD [3] 


inch 


0.983 


18 


11.8 




"uds tag 


0.983 


9 


10.7 




"c tag" 


0.983 


9 


9.6 




"b tag" 


0.983 


9 


9.3 


ALEPH [4] 


inch 


0.97 


13 


11.6 


DELPHI [5 


inch 


1.0 


12 


3.9 




"wds tag 


f 1.0 


12 


0.7 




"6 tag" 


1.0 


12 


9.0 


OPAL [6] 


"m tag' 


1.10 


5 


7.8 




"d tag' 


1.10 


5 


12.8 




"s tag' 


1.10 


5 


5.4 




"c tag" 


1.10 


5 


5.0 




"6 tag" 


1.10 


5 


5.7 


STAR [10] 


P 


0.95 


14 


42.9 




P 


0.95 


14 


32.4 


TOTAL: 






159 


186.1 



of processes with hadrons in the final-state with remark- 
able precision, producing at the same time accurate and 
universal sets of fragmentation functions. The availabil- 
ity of these crucial non-perturbative inputs now opens a 
wide range of possibilities for detailed studies of the nu- 
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TABLE VIII: Same as in Tab. |lV]but now at LO accuracy. 



experiment 


data 


rel. norm. 


data points 






type 


in fit 


fitted 




TPC [2] 


incl. 


1.041 


17 


7.7 


SLD [3] 


incl. 


1.014 


21 


63.0 


ALEPH [4] 


incl. 


1.013 


27 


46.4 


DELPHI [5] 


incl. 


1.0 


12 


12.2 




"uds tag" 


1.0 


12 


22.3 




"6 tag" 


1.0 


12 


16.6 


TASSO [llj 


incl. (44 GeV) 


1.0 


14 


22.3 




incl. (35 GeV) 


1.0 


14 


56.3 


OPAL \12] 


incl. 


1.0 


12 


28.5 




"uds tag" 


1.0 


12 


32.9 




"c tag" 


1.0 


12 


18.1 




"6 tag" 


1.0 


12 


8.6 


ALEPH [4] 


long. incl. 


1.0 


11 


11.5 


OPAL [131 


long. incl. 


1.0 


12 


8.1 


DELPHI [14] 


long. incl. 


1.0 


12 


7.2 




long, "uds tag" 


1.0 


12 


28.0 




long. "6 tag" 


1.0 


12 


5.0 


EMC [18] 


h+ 


1.0 


98 


183.0 




}r 


1.0 


99 


289.1 


CDF [151 


630 GeV 


1.1 


16 


222.4 




1.8 TeV 


1.1 


37 


515.5 


UAl [16| 


200 GeV 


0.988 


31 


84.5 




500 GeV 


0.988 


32 


129.4 




630 GeV 


0.988 


41 


48.3 




900 GeV 


0.988 


44 


391.7 


UA2 [17] 


540 GeV 


0.988 


27 


135.3 


TOTAL: 






661 


2393.9 



cleon structure at NLO accuracy through single-inclusive 



hadroproduction processes, encompassing the study of 
the spin and flavor structure with polarized SIDIS mea- 
surements [26], in proton-proton collisions ^2?], and in 
photoproduction [28|, the modification of parton densi- 
ties by nuclear effects 12911 , or studies of high pT electro- 
production processes [30| . 
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V. APPENDIX: LO RESULTS 

For completeness, we have also performed global anal- 
yses of the same sets of data given in Tables [TT] and IIVI 
where now all observables, as, and the scale evolution of 
the fragmentation functions are computed at LO accu- 
racy. We use the same parametrization ^ and fitting 
procedure as in the NLO case, outlined in Sec. II. The 
parameters of the optimum LO fits for Df and are 
given in Tables |V] and IVII while the values obtained 
for the individual sets of data are compiled in Tables IVIII 
and lVIIll respectively. 

It is important to notice that the total of the LO 
fits is significantly worse than at NLO accuracy, in par- 
ticular, in case of the inclusive charged hadron data. Our 
LO sets should be used only for rough estimates of ob- 
servables where NLO corrections are not yet available, or 
in event generators based on matrix elements at LO ac- 
curacy. Because of the limited usefulness of the LO sets, 
we refrain from going into any further details here. 
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